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Abstract 

The Hawking radiation forms the essential basis of the black hole thermodynamics. The black hole thermodynamics 
denotes a nice correspondence between black hole kinematics and the laws of ordinary thermodynamics, but has been 
so far considered only in an asymptotically flat case. Does such the correspondence rely strongly on the feature of 
the gravity vanishing at the infinity? In order to resolve this question, it should be considered for the first to extend 
the Hawking radiation to a case with a dynamical boundary condition like an expanding universe. Therefore the 
Hawking radiation in an expanding universe is discussed in this paper. As a concrete model of a black hole in an 
expanding universe, we use the swiss cheese universe which is the spacetime including a Schwarzschild black hole in 
the Friedmann-Robertson- Walker universe. Further for simplicity, our calculation is performed in two dimension. The 
resultant spectrum of the Hawking radiation measured by a comoving observer is generally different from a thermal 
one. We find that the qualitative behavior of the non-thermal spectrum is of dumping oscillation as a function of the 
frequency measured by the observer, and that the intensity of the Hawking radiation is enhanced by the presence of 
a cosmological expansion. It is appropriate to say that a black hole with an asymptotically flat boundary condition 
stays in a lowest energy thermal equilibrium state, and that, once a black hole is put into an expanding universe, it is 
excited to a non- equilibrium state and emits its mass energy with stronger intensity than a thermal one. 
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So far various properties of black hole spacetimes have been revealed. Among so many well-known properties of the 
jblack hole, one of the most remarkable ones is the black hole thermodynamics. This tells us a nice correspondence 
• »~j 'between the classical properties of a black hole spacetime and the laws of ordinary thermodynamics || . The black hole 
KnJ thermodynamics is one of the most important stages for understanding the nature of a strong gravity. However this 
^ .nice correspondence has been considered only in the asymptotically flat case. Does such the correspondence rely strongly 
on the feature of the gravity vanishing at the infinity? It has not been clarified if such the correspondence is true of 
the case with a dynamical boundary condition. [|] My general purpose is to extend the black hole thermodynamics to a 
dynamical situation, and hopefully find a non- gravitational system which corresponds to a black hole with a dynamical 
boundary condition like an expanding universe. By the way, the Hawking radiation plays the essential role in the black 
hole thermodynamics to determine precisely the temperature of a black hole with asymptotically flatness ||. Hence, as 
the first step to attack my purpose, this paper is especially designed for the Hawking radiation in an expanding universe. 

There is a serious problem of defining a black hole. A black hole is defined as the spacetime region excluded from 
the causal past of the future null infinity |Ij . The future null infinity for the asymptotically flat case is well known, 
therefore it is rather clear if an asymptotically flat spacetime includes a black hole. On the other hand, however, under a 
dynamical boundary condition, it is generally quite difficult to extract the asymptotic structure of the spacetime directly 
from the Einstein equation. The black hole cannot be defined unless the asymptotic structure is known. In order to 
avoid the problem of defining the black hole in an expanding universe, we make use of the swiss cheese universe. 

The swiss cheese universe (SCU) is the spacetime including a Schwarzschild black hole in a Friedmann-Robertson- 
Walker (FRW) universe. It is not the solution obtained by solving directly the Einstein equation, but the solution 
constructed by connecting the Schwarzschild spacetime with the FRW universe at a given spherically symmetric timelike 
hypersurface E by the Israel junction condition with requiring no surface energy density confined on E || Q. The 
spatial section of E which is a two dimensional sphere is expanding as seen from the black hole side with respect to the 
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Schwarzschild coordinate. So once the two sphere on E was placed outside the gravitational radius R g = 2GM at a 
given time, the event horizon should last to exist after that time. That is, the existence of a black hole is guaranteed in 
the SCU. 

Another technical problem has still been left. Even if with the asymptotic flatness, the curvature scattering makes 
it very difficult to obtain a complete form of the Hawking spectrum in four dimensional case. Indeed, the curvature 
scattering is ignored in the statement that an asymptotically flat black hole is in a thermal equilibrium state. So we also 
want to ignore the curvature scattering for calculating the Hawking radiation in the SCU. By the way, the curvature 
scattering for a minimal coupling massless scalar field vanishes on any two dimensional spacetimes due to the conformal 
flatness of the spacetime. Therefore in this paper, we treat the two dimensional SCU for simplicity, and introduce the 
minimal coupling massless scalar field for calculating the Hawking radiation. 

In section ^J, we briefly explain the two dimensional "eternal" swiss cheese universe and show the strategy to obtain 
the Hawking radiation on two dimensional SCU. The behavior of a minimal coupling massless scalar field on the eternal 
SCU is analyzed in section |3|. Section |^ is devoted to computing the Hawking radiation. Summary and discussions are 
given in section ^|. 

The unit throughout this paper isc = h = ks = 1. The metric signature is (—,+,+,+) for four dimension and 
similarly for two dimension. 



2 "Eternal" swiss cheese universe and the strategy 
2.1 Two dimensional "eternal" SCU 

The two dimensional SCU should be constructed by connecting the Schwarzschild spacetime with the FRW universe at a 
given timelike surface (curve) E. The way of junction is very similar to the four dimensional case explained in appendix 
|Aj. Hereafter as the terminology, let the word "BH-side" denote the spacetime region inside E where the metric is given 
by (|]) or (||), and "FRW-side" for the region outside E where the metric is (^) or ((?]). If a quantity Q is measured in 
the BH-side, we express it as Qbh and similarly Qp in the FRW-side. 
The metric of the BH-side is given by 

ds 2 BH = — (l — dT 2 +(l- dR 2 , in TR-system (1) 



R J V R , 
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= e - R/R <> [~dW 2 + dX 2 ] , in WX-system , (2) 

where R g = 2GM, G is the gravitational constant, and the meaning of M is given at the end of this paragraph. 
The terms "TR-system" and "WX-system" mean respectively the Schwarzschild coordinate and the Kruskal-Szekeres 
coordinate in the BH-side. Because the maximally extended Schwarzschild spacetime has two exterior regions of a black 
hole, there can also be two exterior regions in the BH-side of the SCU. We distinguish them by calling "L-region" and 
"R-region" . The coordinate transformation between the TR-system and the WX-system is given by 

W-X = ^fexp[-K(T-R*)} , W + X = iexpKT + i?*)] , (3) 

where the upper signature is for the R-region and the lower for the L-region, R* = R + R g \n{R/R g — 1) and k = l/2R g . 
We set the Killing vector (;bh in the WX-system as 

(,bh = dw , all over the BH-side , (4) 

where a tilde is added to the quantity measured in the WX-system. Adopting the direction of this Killing vector as the 
standard of the future direction, the Killing vector ^bh in the TR-system is given by 

A I 8t , in the R-region 

£ - BH = 1 a i -i j . (5) 

I —Ot , m the L-region . 

The normalization of the vectors (Q) and (g) is determined at the asymptotically flat region of the full Schwarzschild 
spacetime which is not present in the SCU. However once we accept them as the timelike Killing vector, k can be 
interpreted as the surface gravity of the event horizon, and further M becomes exactly equal to the Komar mass |jj 
evaluated on any spatial section of E. Such a Komar mass Q defined using our Killing vector (|4|) or (||) is invariant 
under the deformation of the spatial section of E. See the end of appendix |A] for more detailed explanation of k and M. 
The metric of the FRW-side in the comoving coordinate is 

ds 2 F = -dt 2 + a(t) 2 -^—, (6) 

_1_ fZT 
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Future Event Horizon 



Past Event Horizon 



FRW-side in L-region 




Junction Surface 



FRW-side in R-region 



Future Direction of t and T 
along Junction surface 



Figure 1: Two dimensional eternal swiss cheese universe. The junction surface is set outside the gravitational radius so that a 
Black hole exists eternally in the BH-side. 



where k — ±1, is the spatial curvature, and a(t) is the scale factor. The comoving spatial coordinate in two dimensional 
case can range over infinity — oo < r < oo for open and flat cases k = —1,0. There is coordinate singularities at r = ±1 
for closed case k — 1, but it can be eliminated by the transformation from the comoving coordinate (i, r) to the conformal 
coordinate r*) which is given by dr\ = dt/a(t) and r = sinr*, r* , or sinhr* for k = 1, or — 1 respectively. Then the 
metric becomes 

ds% = a^f (-dry 2 + dr* 2 ) , (7) 

where — oo < r* < oo. Hereafter as the terminology we call the conformal coordinate the Cfl-system. The comoving 
coordinate is equivalent to the Cfl-system in the sense that there exists one comoving coordinate point (t, r) for each 
point of the Cfl-system (r],r*), where this correspondence is the one-to-one mapping for open and flat cases and the 
onto-mapping for closed case. 

The two dimensional FRW metric does not have any timelike Killing vector, but has a conformal timelike Killing 
vector £p which satisfies £p + oc g^, 

£ F _ f 9 n ' in tne R - re g ion ^ 
1 —dri , in the L-region . 

In order to treat the two dimensional SCU as a simplification of four dimensional case, we require the same kinematics 
of the junction surface and the same relation between the temporal coordinates t and T as four dimensional case shown 
in the appendix [A|. That is, the location of S is given by 

| R(t) = a(t) r a , in the BH-side, ^ 
1 r = r (= const.) , in the FRW-side, 



and the temporal coordinates are related by 

dT(t) _ 177 (1Q) 



dt 1 - R g /R(t) 

The unit normal vector to S is given by the same form as the four dimensional case (^) and ( |(34| ) 

rip — ( , — t-t ) in Cfl-system, (11) 



o(t). 

R(t) ^ Rg 



BH 



1 - R g /R{t) ' V R(t) 



t{t) in TR-system, (12) 



where Q = dQ/dt. Note that the equation ( |To| ) and the normalization of (p"2|), rig H = 1, gives the Friedmann equation 
of dust matter (|66[). However in order to extract the essence of the Hawking radiation in an expanding universe in a 
simple way, we dear to assume that the scale factor satisfies a(t) ro > R g for — oo < t < oo. That is, we consider the 
two dimensional SCU in which a black hole can exist eternally. The conformal diagram of such the "eternal" SCU is 
shown at the figure [l]. With this assumption, a black hole can be defined without respect to the spatial curvature of 
the FRW-side, and the event horizon bifurcates. Hereafter throughout this paper, we adopt the eternal SCU as the 
background spacetime. The validity of the eternal SCU will be discussed at the end of section ||. 
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2.2 Strategy for the Hawking radiation 



2.2.1 Choice of the Vacuum state 

As discussed in section |l], we introduce a minimal coupling massless scalar field $ and neglect any back reaction of <!> to 
the background eternal SCU. The purpose of this paper is to compute the Hawking radiation received by a comoving 
observer in the FRW-side. Then, because of the equivalence between the comoving coordinate and the Cfl-system, the 
expectation value of the number operator of the quantized scalar field in the Cfl-system attracts our interest. 

Here we have to specify on which vacuum state the expectation value should be calculated. For the asymptotically 
flat black hole spacetime, there are two candidates for such the vacuum state: the Hartle-Hawking state 11 and the 
Unruh state The former describes that a black hole is surrounded by a thermal cavity of temperature k/2tt and 
the whole system is in equilibrium. On the other hand, the latter corresponds to the situation that the environment 
surrounding the black hole is not a thermal cavity but an empty space, that is, the black hole is evaporating due to the 
Hawking radiation without absorbing any energy from the empty environment. By the way, it is not a priori known 
whether the black hole in the SCU emits a thermal radiation or not. Therefore it seems appropriate for the SCU not to 
put the black hole in a thermal cavity, but to put in an empty environment. Then we search what kind of spectrum is 
radiated out of the black hole in the SCU. From above considerations, we adopt the Unruh state as the vacuum state 
on which the Hawking radiation is computed. Further due to the presence of the junction surface, we use the following 
strategy of computing the Hawking radiation. 



2.2.2 First step: construction of the mode functions 

The vacuum state is constructed by quantizing the scalar field <i>. In order to quantize <£>, we have to obtain the positive 
frequency mode functions. On the SCU, there are two categories of mode functions: one of them is a collection of 
the mode solutions which is monochromatic in the BH-sidc but a superposition of some monochromatic modes in the 



FRW-side because of the junction conditions of <&, (|68| ) and (|69| ) given in appendix |B|. Another category is of the 
monochromatic mode in the FRW-side but a superposition in the BH-side. And note that a mode function generally 
takes a different functional form according to the coordinate system in which the Klein-Gordon equation is solved. Then 
we search the following three modes; 

1. Monochromatic mode in BH-side with respect to the WX-system, which we call the WX-mode 

2. Monochromatic mode in BH-side with respect to the TR-system, which we call the TR-mode 

3. Monochromatic mode in FRW-side with respect to the Cfl-system, which we call the Cfl-mode 

Because our background is two dimensional and $ is massless and has no coupling to gravity, the curvature scattering 
of $ never occur. Therefore what we have to take care is the possibility of the reflection of $ at S. Further as seen in 
the next section, it will be revealed that no reflection at £ arises due to the special property of our junction conditions, 
and we can accomplish completely the construction of these modes. Then we define the positive frequency WX-mode, 
TR-mode and Cfl-mode, with respect to the (conformal) timelike Killing vectors (Q), (^) and (||). That is, the positive 
frequency mode should satisfy = — iuQ, where uj(> 0) is the frequency. 



2.2.3 Second step: quantization of the scalar field 

We quantize the minimal coupling scalar field $ in the WX-system, the TR-system and the Cfl-system using the mode 
functions obtained in the first step. Then the vacuum states of each mode functions are defined. 

The observer is comoving in the FRW-side, so the vacuum state appropriate for describing the empty environment 
surrounding the black hole is that of the Cfl-mode. Therefore the Unruh state for the SCU should be composed of two 
vacuum states: one of them is the vacuum of the WX-mode \0} w on the past event horizon and another is that of the 
Cfl-mode |0) ?J on the past null infinity Q 

By the way consider if the background is four dimensional, then it is very difficult to solve the Klein-Gordon equation 
on the whole spacetime due to the curvature scattering. However near the event horizon and the null infinity, the solution 
can be found. Therefore for four dimensional case, we cannot help using the Unruh state as our vacuum state on which 
the Hawking radiation is calculated. However on the two dimensional eternal SCU, the exact form of the mode function 
can be constructed on the whole spacetime in the first step. So due to the absence of the curvature scattering and the 
reflection at S, the out-going modes radiated from the past event horizon are received by the comoving observer without 
loss of the energy of the out-going flux, and the in-going mode from the past null infinity never reach the observer. That 
is, the vacuum 10)^ on the past null infinity has no effect on the Hawking radiation received by the comoving observer. 
Consequently it is enough for us to prepare the vacuum state \0} w at one Cauchy surface as the vacuum state on which 
the Hawking radiation is calculated. Hence in the following sections, we do not take care about the vacuum (0)^ in the 

2 For asymptotically flat case, the observer rests on _R = constant, so the vacuum of the TR-mode is chosen on the past null infinity. 
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0) at Cauchy Surface Comoving Observer 



0) ,„ at Past Event Horizon 




Reflection at £ ? 




Reflection at£ 
and Curvature Scattering 



0) at Past Null Infinity 



Figure 2: Right figure shows the vacuum state for four dimensional case. Right one is of two dimensional case. In section ^ it 
will be found that no reflection takes place for two dimensional case, and it is enough to prepare the vacuum \0) w on a Cauchy 
surface for computing the Hawking radiation in the two dimensional SCU. 

Unruh state, and aim to calculate just the Bogoljubov transformation from the WX-mode to the Cfl-mode. This can be 
understood as the particle of the Cfl-mode created in the vacuum state of the WX-mode. 

To proceed this calculation, we make use of another intermediate transformation, that is, we perform the successive 
Bogoljubov transformation from the WX-mode to the TR-mode, and from the TR-mode to the Cfl-mode. Using this 
successive transformation, the difference between the asymptotically flat case and our SCU case should be clarified. The 
first transformation from the WX-mode to the TR-mode gives the thermal Hawking spectrum since the bifurcate event 
horizon of the eternal SCU lets us to make the same discussion given in the reference 0. Then the deviation from the 
asymptotically flat case should arise from the second transformation from the TR-mode to the Cfl-mode. 

3 Massless scalar field in the eternal SCU 

3.1 Classical mode function 
3.1.1 Preparations 

We seek the positive frequency TR-mode, WX-mode and Cfl-mode on the SCU. However a positive frequency mode is 
out-going for positive wave number, and in-going for negative wave number. Therefore we look for the out-going and 
in- going modes. Those modes can be constructed by connecting the mode function obtained on the full Schwarzschild 
spacetime with that obtained on the full FRW spacetime at E. Because $ is a minimal coupling massless field, it is 
natural to require that no potential of $ is confined at S, that is, the field equation of the form □$ = holds even at 
E. This requirement determines the junction conditions of $ at E, 

= $FBW| E , dn^BH^ = <9„$Fw| S , (13) 



where d n = n^d^ and is the unit normal to E, ( |ll| ) or (|12|). The detailed derivation of these conditions is given in 
the appendix [^. 

The mode functions satisfying □<£> = on the full Schwarzschild spacetime are, 

,(±) -•_ ( e -in(T-R*) out _going . 

f„ff = — ; x < D *s m IR-system, (14) 

BH ' n |e- i0 ( T +«) , in-going ' 

7(±) 1 ( e -ia(W-X) ; out-going . 



where the upper suffix "+" denotes the out-going mode and the lower "— " is for the in-going one, and the normalization 
factor is determined by an ordinary inner product of the full Schwarzschild spacetime. The mode functions on the full 
FRW spacetime are 



1 J e * u w r ' , out-going 



v/SttH I e- lw W+ r > , in-going 



in Cfl-system , (16) 



where the upper suffix "1" denotes the out-going mode and the lower "2" is for the in-going one, and the normalization 
factor is determined by an ordinary inner product. Even if we solve □$ = using the comoving coordinate, the same 
mode function is obtained. For later convenience, we introduce the notations: 

/ F)W (^) = e-^ and h^\r*) ee —L— * . (17) 

V 47r M 
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Note that, because of the assumption of the eternal SCU, the temporal part fF,u{rj) can span a complete orthogonal set 
with respect to 77 G (—00, +00) and lo S (—00, +00) as 



d V /f> (v) /f>' fa) = 2tt <5(u> - a/) 



/*>(»/) =27r«(>7-i/) 



(18) 



Referring to the ordinary definition of the inner products on the full Schwarzschild and FRW spacetimes, we define 
that on the SCU, 



K*(>?) 



dR*ipld p ip 2 + i / dr*i>ld p ip- 



(19) 



where ip denotes the solution of □$ = on the SCU, d p = (l/I^Q^c^, £ denotes the (conforms!) Killing vectors (||) and 
(||) or (||), and is the coordinate value of r* at the junction surface S. This inner product is defined on the spatial 
surface of 77 = constant in the SCU. 

3.1.2 Out-going and in-going TR-modes and WX-modes 
TR-mode 

As mentioned at the subsection 2.2. 2| , the out-going and in-going TR-modes are monochromatic in the BH-side but a 
superposition of some monochromatic modes in the FRW-side. Then we can express the out-going TR-mode V'bh n ana - 
the in-going one ip BH \ q as 



tBH,Q 



<Pbh,ci 

^BH,Q 

f°° did 
J —00 



.4 



(i) ^(i) 



A 



(2) ,(2) 
Qu Vf,u 



(i) ^(i) 



5, 



(2) . (2) 



, in BH-side 
, in FRW-side , 

. in BH-side 
, in FRW-side , 



(20) 
(21) 



where the junction coefficients A' 1 ' 2 ) and Z?' 1,2 ) are determined by the junction conditions (|T3|), which become 



dio 



A ( nl Mv) ft£i(r ) + ^ /*>(»?) ftgLfo) 



> {1) 
'f. 



(2) 



,(2) 



1 /*°° 



(1) 



(2) 



Mv) "^T M + A nl Ml) ~£r(ro) 



dr* 



(22) 



and the similar equations for 2?"' ' holds, where (T(r)), R(r))) is the coordinate point on £ in the TR-system, and note 
that (dhF / dr*)\ r — ro = constant. 

Here in this subsection, let us calculate to obtain A^ 1,2 \ then the similar calculation gives B^' 2 \ Using the com- 
pleteness (|l8|), we obtain the integral representations 



A 



(i) 



A 



(2) 



sgn(w) 
sgn(w) 



e iujr* 



4?rv1^J| 



[i a (si,u) - j A (n,w)] 
[i A (n,Lo) + j A (n,io)} 



(23) 



where sgn(aj) = w/|u;|. Here I A and J A are given by 



I A {n,u)=w drje ^r,-in[T{r,)-R*(n)] ^ j A (n,u;) = f2 



dr)P A (r]) e 



iwr)-in[T(ri)-R* (77)] 



(24) 



where 



^(r?) 



0/(77) r 



d?7 ' 



and a' = da/drj. Once the scale factor is specified, the junction coefficients can be determined in principle. One may 
think that the out-going TR-mode is a superposition of the out-going and in-going modes in the FRW-side, and so is 
the in-going TR-mode. However it is not the case. To show it, it is essential to notice that 



dr) 



e , M ,-,!l[T(,)-K*(,)] = i r^ + QP^fe)] e «-»l-itolT(ri)-R*(r))] 



(25) 



G 



which gives 

I a + Ja = -% e l 
Then we obtain the out-going TR-mode in the FRW-side to be 



_ _• iuj V -iO.[T{ri)-R*(ri)} 



77—OC 
T}= — OO 



(+) 

BH.n 



in FRW-side 



-in[T(fj)-R*(n)] 



(4tt) 3 / 2 VN 



77—OO 
77— — GO 



(26) 



(27) 



where <Ae (1/x) exp(±?x) = ±i7r is used. As discussed in the appendix |c| the second term should be zero, and we 
can find I a — — J a due to the equations (p6|). 

The same way of calculation done so far holds for the other junction coefficients B^' 2 \ Then the out-going and 
in-going TR-modes are obtained to be 



rBH,n 
tbh,q. 



, in BH-side 
, in FRW-side , 

, in BH-side 
, in FRW-side , 



where A^> and B^ are summarized as 



27rV|^w 



2tt 



<i?7 e 



<wr7-ifi[T(77)-fl*(?7)] 



(28) 
(29) 

(30) 



and 



R (2) _ e 



2tt 



drj e 



iwri-iCl[T(n)+R*(ri)] 



These results mean that the junction surface does not reflect the scalar field, therefore the out-going TR-mode is a 
superposition of the out-going modes in the FRW-side without including any in-going mode, and the in-going TR-mode 
consists of the in-going modes in the FRW-side. Further we can check by some routine computations that the inner 
products of these modes defined by ( |l9| ) result in 



(±) 

BH,W 



sgn(fi) S(Q - SV) , , ^ n ,) = . 



(31) 



This means that the out-going and in-going TR-modes given by (|28| ) and (|29|) span a complete orthonormal basis. 

There are two reasons for the disappearance of reflection at E: (i) the junction condition for the scalar field with no 
surface potential at E, and (ii) the conformal flatness of the background spacetime. To understand the first reason, let 
us consider the case of F 7^ for the junction condition ( |69| ) in the appendix |b[ In this case we can follow the same 
way of calculation to obtain a similar result to 



(i) 

(2) 



/oo 
drj a fp u F 
-00 

/oo 
d/qa fp u F . 
-00 



The second term of A^> does not vanish, and the reflection should take place at E. On the other hand even if F = 0, 
when the background spacetime is not conformal flat, cannot vanish due to the curvature scattering. 



WX-mode 

The out-going and in-going WX-modes 5 can be obtained by a similar way to the TR-mode. The resultant 
WX-mode is not reflected at E too, 



(+) 

BHS1 



$<+> - 

BH,n 



, in BH-side 



f 00 dw A { } ] <$;], , in FRW-side . 



(-) 

BH,n 



BH,n 



r dcuB {2) ^ 



, in BH-side 
, in FRW-side . 



(32) 
(33) 
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The junction coefficients A^ and B^ are given by 




ijW = e <-[ J I dr] e lujr > ' " ! 1 ' '' 1 " N 1 " ' 

where (W(rf), X(r])) is the coordinate of a point on £ measured in the WX-system. These WX-modes $^1 ^ also span 
a complete orthonormal basis. 

3.1.3 Out-going and in-going Cfl-modes 

(12) 

The out-going and in-going Cfl-modes ip F ^ are monochromatic in the FRW-side and a superposition of some monochro- 
matic modes in the BH-side. Because the reflection of $ at £ disappears as discussed in the previous subsection, these 
modes can be expressed as 



(+) A+) 




, in BH-side 



V^l = J- T ai ^n^BH,n (34) 
1 J! , in FRW-side, 

, in BH-side 
, in FRW-side . 



= v,( 2 r (35) 



where we use the TR-mode in the BH-side for later use. The junction coefficients and D' _ ) can be determined 

in principle by the junction conditions of the scalar field. However because the junction surface does not rest in the 
TR-system R(rj) = a(r/) ro, it is technically difficult to follow the same calculation done in the previous subsection. In 
other words, the completeness of the temporal part, e _lOT , cannot extract the simple integral representation of the 
junctio n co efficients C 1 -^ and like (p3|). However the analytic representations of them will be obtained in next 

section 3.2 with making use of the Bogoljubov transformation. Therefore we assume temporarily that these Cfl-modes 
span a complete orthonormal basis. After obtaining the integral representations of C^ + ' and D^~\ we will find that this 
assumption is consistent. 



We can easily find some relations between A (1 \ and C^ + \ \ By substituting ( |28| ) into (34) and other 
similar operations, we obtain 

/oo />oo 

duA^c™ = *(n-n') / dnig ) cw=j( w - w ') (36) 

-oo J —oo 

dwB£lrf-l = 6(n-n') / dnB^rf-J, =6{u>-tJ) . (37) 



3.2 Quantization 

3.2.1 Vacuum states of WX- mode, TR-mode and Cfi-mode 
Positive frequency modes 

For canonical quantization of the scalar field <&, we need the positive frequency mode. Hereafter as the notation, we 
denotes the positive frequency TR-mode, WX-mode and Cfl-mode respectively as ^bh, ^bh and ^p- These modes 
can be obtained from the out-going and in-going modes as 



ib for K > 

^bh,k = { f"y n=K \ T ^ „ (Positive Frequency TR-mode) 



,/+) 

BH 

^BH.n=-K ,fovK<0, 
where f2 is the frequency, K is the wave number and f2 = \K\. Similarly we find 



and 



fV4 + 2o & ,fotK>0 
^BHK={ f H ,' ~ (Positive Frequency WX-mode) 



7^= fc ' (Positive Frequency Cfl-mode) 

^ F L=-k ,forfc<0, 
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where = \K\ and u> = \k\. For example, the positive frequency TR-mode in the BH-side is expressed as 



v f bh.k in BH-side : 



-HQT-KR*) 



1 < f e -«r(T-H*) , for X > 
X 1 e+ iif ( T+iJ *) , for K < , 



and this satisfies the definition of the positive frequency mode C^ BH ^ bh = —iQ^bh where > 0. 

We should note that the mode functions ^>bh and >J> are defined only in the R- region or in the L-region of the 
maximally extended eternal SCU. Therefore we define the other mode functions as M 



R^ BH.K — 




, in R-region 
, in L-region 



BH.K — 



, in R-region 

^bh.k , in L-region, 



and 



R^F.k 




, in R-region 
, in L-region 



F.k 




, in R-region 
, in L-region . 



All of these modes are of the positive frequency, and can be defined on a Cauchy surface connecting the spatial infinities 
of R-region and L-region via the bifurcation point of the event horizon. We take the complex conjugate in the definition 
of l^bh.k and F.k, since the future temporal direction in the L-region with respect to T and rj is inverted in 
comparison with that in the R-region. Further it is useful to define the following mode functions [9j, 



1 



II 

BH.K 



^sin^Trn/K)! 1 ^ 
1 

|2sinh(7rlV0| 1/2 



e R^BH.K + e L^BH.K 

e R™BHK + e L^BH.K 



(38) 



These mode functions share the same analyticity as the mode S&bh- That is, ^ B ^i and ^ bh are analytic and bounded 
on the lower-half-plane of the complex U{— W ± X) plane. 

Due to the completeness of the out-going and in-going TR- modes ( |3l| ) , the positive frequency TR- modes satisfy the 
complete orthonormal relations. So is the positive frequency WX-modes. For the positive frequency Cfl-modes, the 
assumption mentioned at the subsection 3.1.3 makes the positive frequency modes satisfy the complete orthonormal 
relations. Hence we can carry out the canonical quantization of $ using each mode function. 

Quantization with WX-mode 

The scalar field <£> is quantized using the WX-mode as 



$ = 



dK 



~ a K ^ BH.K 



K BH.K 



Further because '5 bh and '5 bh share the same analyticity, we can find another representation, 



dK 



a K H H ,K + 4 f *BH,K + a K *BH,K 



BH.K 



Here a^> and a^ 1 are the annihilation operators, and their Hermitian conjugates are the creation operators. Due to 



the same analyticity of ^ B li as ^flfli they also share the same vacuum state \§) w defined by 



a K \0) w = |0> 



ir 



|Q) 



ir 



for all K and K . 



Quantization with TR-mode 

The quantization of <J> using the TR-mode is represented as 



$ = 



dK 



R&K R^ BH.K 



BH.K 



+ L b 



L»K L 1 ^ 'BH.K i" L°x L * BH.K 



+ L bl 



where nbx and i}>k are the annihilation operators and r$ k and i$ K are the creation ones. The vacuum state of the 
TR-mode |0) T is defined by 



Rb K |0) T = L b K |0) T = for all K . 
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Quantization with Cfl-mode 

The quantization with the Cfl-mode is similar to that of the TR-mode, 

/oo 
dK R C k R V F ,k + R c{ R ^ F ,k + L C k L ^BH,K + Lc\ L^* F J 
-oo 

xc the 

Cfl-mode |0)^ is defined by 



where R c k and ^c k are the annihilation operators and R c\ and lc\ are the creation ones. The vacuum state of the 



RC k \0) n = L c k \0) v = for all fc . 

3.2.2 Bogoljubov transformation 

The Bogoljubov transformation is the change of positive frequency basis. Therefore the equation ( |38| ) gives the Bo- 
goljubov transformation between the WX-mode and the TR-mode. In this subsection, we derive the Bogoljubov trans- 
formation between the TR-mode and the Cfl-mode. It is important to notice the implication of the relations ( |36| ) and 
( |37| ) that, for the out-going modes, 

/>oo 

d«4!Mln = / MC$L ) <t>£l a = $l inBH-side 

-OO 
/• X 



dJ dfl C { +J An*' 4>%, = in FRW-side , 

and a similar one for the in-going mode. This denotes that the junction coefficients B^ 2 \ C*W and £)( — ) give the 
change of out-going and in-going basis as 

/•OO P oo 



J — oo J — oo 

and 



/OO pOO 
-oo J — oo 



n ■ 



Consequently, by definition of ^ bh, ^ f, ^bh ano - '0f ' 2 ' > i we obtain the Bogoljubov transformation between the positive 
frequency modes $m and ^ f as 



/oo 
dk [a Kk t> F ,k + &Kk ^* F ,k) 
-OO 

/OO 
dK [X Kk ^BH,K + ilKk ^*BH,k] 
-oo 



where the Bogoljubov coefficients are given as 



x Kk 

otKk = { , for K > and k < 0, or K < and fc > 



A^j. , for K > 0, fc > 



B-K-k ,for ^ < 0, fc < 0, 
, for K > 0, fc > 



(39) 



3 K k = {0 , for K > and fc < 0, or K < and fc > 



i?i 2 L, , for K < 0, k < , 



and 



'cj+ } , for if >0, fc > 

S " , for if > and k < 0, or K < and fc > 

£>^t_ fc ,for if < 0, fc < 0, 

Ci^ fc ,forif>0, fc>0 
fJ-Kk = {0 , for K > and fc < 0, or if < and fc > 

L>^ fc ,for if < 0, fc < 0. 
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Further, taking the inverse of the Bogoljubov transformation, we can find the relation between (a, (3) and (A, fi) to 
be \ K k = a* Kk , fi K k = -f3* Kk - Hence the Bogoljubov transformation from (r^bh.k , l^bh,k) to (r^f^, L^F.k) is 
obtained to be 

/oo 
dk [a Kk R^ F ,k+Pi<k R^*F,k] 

(40) 

L^BH.K = 



/oo 
dk [a* Kk L^F,k + P*Kk L*F,fc] , 
-oo 



Finally note that the junction coefficients and are represented analytically using the other junction 

coefficients A 1 - 1 ) and B^> through the Bogoljubov transformations obtained above. Therefore we can find that the 
out-going and in-going Cfl-modes span a complete orthonormal basis. 



4 Hawking radiation in the SCU 



4.1 Hawking radiation 



We compute the Hawking radiation observed by a comoving observer 

4 



As discussed at the section 2.2 



its spectrum is 

given by w (0| Np M \0} w , where Np^ = rc}^ rc u is the number operator of particles of the Cfl-mode, [] and lo(> 0) 
is the energy of the particle which is equivalent to the frequency of the positive frequency Cfl-mode. Further in order 
to treat the situation that the observer receives the out-going radiation from the black hole, we should consider only 
the out-going Cfl-mode, k > 0. Then taking the Bogoljubov transformations (|3^) and ([l0]) successively, the Hawking 
spectrum is obtained to be 



w 



<o| Np iU \o) w = J™ dn e2 J lK _ x D H (Q, 



where Dh is given by 



D H (Q,u) 



A 



(i) 



A 



(i) 
si -u 



2e 



lS!/fi s(4! 



A 



(i) 



(41) 



(42) 



and 5i denotes the real part, A^ is given at (^0|), the relations fl = \K\ = K and lo = \k\ = k are used, and due to the 
relations (|3|) the variable of integration is restricted to K > fl > 0. 

Obviously the Hawking radiation in the SCU is totally different from a thermal spectrum. Q The factor Dh represents 
such a difference, and we call Dh the deviation factor. Note that the particle creation due to the cosmological expansion 
does not take place with a minimal coupling massless scalar field in two dimensional background spacetime jToj] . So 
the deviation factor Dh includes only the effect of the motion of the observer relative to the TR-system, and does not 
include the cosmological particle creation. 

The other notice of the spectrum (^Tj) is that the Bogoljubov transformations used in deriving the Hawking radiation 
are defined all over the spacetime. Especially whole of the temporal information from r\ = — oo to r; = oo are included 
in this result. That is, we should interpret the observer is at a remote future r\ — * +oo. 



4.2 Divergence and normalization of the Hawking spectrum 

The black hole should lose its mass energy M due to the Hawking radiation. However, since the back reaction to the 
background spacetime is ignored in our derivation of the resultant spectrum ( f4l| ) , the mass M remains constant during 
emitting the Hawking radiation. Thus, the spectrum (^) may give us the conclusion that the observer at a remote 
future receives an infinite energy from the black hole via the Hawking radiation. In the next subsection 4.2.1, 



we see 



the divergence of the spectrum mih for the simplest case that the scale factor is constant. Then, the normalization of 



the spectrum (41) is discussed in the following subsection 4.2.2 



4.2.1 Simplest case: constant scale factor 



Consider the case a - 
then (0) gives T(rj) 



a c = constant, which should give the thermal spectrum. From (J9|), R(rf) = a c i*o = constant, 
b c r]+ constant, where b c = a c yl — k r 2 /(l — R g /a c r Q ). That is, T{rf) — R*(r/) = b c i] — R*, 



3 Here we set the observer in the R-region. But even if the observer is in the L-region, the resultant spectrum is the same. 
4 But the case that a = constant should reproduce the thermal spectrum. 
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where R* = constant. Then the junction coefficient given at ( p0| ) is calculated using the representation of the delta 
function dk exp(ikx) = 2ir5(x), 



A 



(i) _ 



S(uj — b c fl) 



(43) 



Since = for Q > and lo > 0, the deviation factor becomes Dh = \A^l\ 2 which gives the spectrum ( [ill ) as, 

b c 



w(0\N Fju \0) w 
This spectrum diverges by the delta function S(0). 



5(0). 



(44) 



4.2.2 Finite time interval normalization 

If the delta function appeared inJ0) is removed, the divergence in the previous resultant spectrum (|4^ ) is normalized. 
Note that the delta function in (|43|) comes from the integral of infinite time interval in the representation (|30|). Therefore 
we make the comoving observer receive the Hawking radiation during a finite time interval, ?/; < r\ < f)f(= rji + rj p ). 
Further we introduce the periodic boundary condition for the scalar field, &(rj + r] p ) — $(77), which changes the frequency 
of the Cfl-mode u> and that of the TR-mode O with the discrete quantities, 



2tt 
lu = — n 
rip 



n = ±1, ±2, ±3, • • • and rj p — r\f — rji 

iV = ±l,±2,±3,--- and T p = T(r) f ) - Tfa) , 



(45) 



where, because the mode solution of zero frequency is trivial $ = constant, the zero frequencies n — and N = are 
neglected. The integral with respect to ui, fi and rj should be replaced as 



J — t 



du) 



E 

n— — oo , ^0 



dn 



E 



drj 



'it 



dj], 



and one of the equation in the completeness ( pl| ) is changed to the form J Vf drj fp u (if) fF, u > (r/) = Vp <^,w' ■ Consequently, 



the junction coefficient A^ can be calculated by the similar way to obtain (|30| 



.4 



(i) _ 



if 



drj e 



iu-q~in[T(r))-R* (r})\ 



(46) 



and A^l — should also hold. Finally the Hawking spectrum (f4l|) is modified to the following form, 



oo _^ 



(47) 



N=l 



where u> and are given by (f45|), and the deviation factor D^f(Sl,w) is given by (|42|). 

Applying this finite time interval normalization to the simplest case of the previous subsection 4.2.1 , the delta 
function in (^) is replaced by the Kronecker delta <5 Wj h c o- Then the spectrum ( ff~4| ) is normalized to the form (f47|), 
I H { m ,rj f -L>) = b c (e 2 ™/ b ^-\)-\ 

4.3 Some models 



In order to understand the properties of the normalized Hawking spectrum in the SCU (47), it is useful to consider some 
cases of the phase factor, T{ij) — R*(ri), in the representation (flq). 

4.3.1 Case: the thermal spectrum 

Here we look for the scale factor a(rj) which makes the spectrum (|47]) be the thermal one. According to the previous 
section 4.2, the phase factor, T(rj) — R*(rj), should be linear of 77, 



T(rj) — R*(rj) = brj + p , b = const > , p = const. 



(48) 
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The spectrum ( |47| ) of this case becomes thermal, Ih{th, rjf, u) — 7j( e 27I W bK — l) 1 , The differential of the phase factor 
(Eg) gives the following, 



da(f]) _ a/1 - kr 2 b ( R g 



<V) - - 1 - -7±r ■ (49) 
dry r r \ a(rj)r J 

Though the simplest case, a = constant, is the trivial solution of this equation, the general solution is a dynamical 
universe, a ^ constant. 

For sufficiently small b, we can set da/di]]^ > at an initial time 770, and generally the scale factor grows exponentially. 
If we set roughly 0,(77) = e an , where a is a constant, then it is rewritten to the form a(t) cx t through the relation 
dt = 0,(77) dr). This is the asymptotic form of the scale factor for the open FRW universe at a remote future t — > +00. 
That is, one may think that the black hole in the open universe emits a thermal Hawking radiation. However the 
situation is not so simple, as discussed in the following two subsections. 

4.3.2 Case: T(rj) - R*{q) cx rj 2 



In this subsection we treat a case that the scale factor gives the phase factor of the junction coefficient (46) as 

T(rj) — R*(rj) = cr/ 2 + bt] + p , c = const > , b,p = const. 
Taking the differential of this equation, we find 



dr) r Q "' r \ 0(77) r 

Generally for the positive initial value 0(770) > at an initial time j]q < —b/2c, this equation makes the scale factor 
0(77) grow exponentially. Therefore by the same discussion given in the previous case, the present case also includes the 
asymptotic form of the open FRW universe at a remote future. Though the rough estimates of the equations (|4^) and 
( p0| ) gives the same conclusion that a(t) cx t, but the difference of them is the detailed behavior of a{t). For example, 
because c > 0, the velocity of the expansion da/drj of the present case is smaller than that of the previous case for 77 > 0. 

The spectrum ( |47| ) of the present case is expected not to be the thermal one, since the phase factor, T{rj) — R*(i]), 
is different from that of the previous case. The normalized junction coefficient (46) is rewritten using the change of the 
conformal time, 77 = Vcfi e l7r / 4 [77 — (uj — bQ)/2ctt], 



4^2 = -77 \/— / dve^, (51) 



e 



~*Pc \,,,\ fVt 



Vp 



n 



where we set £2(= 2ttN/T p ) > 0, that is, N = 1, 2, 3, • • • , and p c = ojr„ - i(ui - bil) 2 /4:cft + pfl + 7r/4. Because the 
integrand has no pole on the complex 77-plane, the integral in this representation of becomes 

fjf r r Mm) 

= -/-/-/ . (52) 

where the integration paths are shown at the figure [| Note that the path L2 is given by 77 = $i(fji) + iy for < y < $s(fji) 
where 3 denotes the imaginary part, then we find 

dfje-v 2 =i e -^i? ( (h) dy e y 2 - 2m ^)y — >0 as 1^1 — > 00 , 



where 5R(?7i) 2 = (cO/2) [77, — (u> — bil) /2cf2] 2 . Similarly it is found that J L — > as \r)i\ — > 00. Therefore we obtain the 
approximate form of A^> for sufficiently large rji, 

42^— o V M / df i e ~ f ' =-7r\/- [Erfc^-Erfc^))], 
■q p tt V c J^ m) ripll V c 

where Erfc(z) = due~ u2 is the Gauss's error function. This gives the deviation factor fl42| ) for large rji, 

D H (Q,Lu) ~ — ^ [AE(fi, W ) 2 + e 2lr °/«AE(f7,- W ) 2 



+2e 7r ° /K cos f b2 ° 2 ^ CJ " - 2p Q - ~ ) AEfO, w) AEfO, - w) , (53) 
4c s 2 
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Im n 



4 \ 



if 

Li 
- Re tT 



Im u 



Reu 



Figure 3: The left figure shows the integration path for (B2) in the 77- plane, and the right one for (|56[) in the it- plane 



where AE(fi,w) = Erfc(5R(^i))-Erfc(5R(%)), = 2irN/T p , uj = 2nn/r] p , and7V ; n = 1,2,3, • ••. Note that 5ft (77) depends 
on f2 and w by definition of 77 given just before (pi]). 

Here the following relation is important; Erfc(a;) = (1/2^/2;) e _2; I 2 W_i/ 4 ,i/ 4 (x 2 ), where W M)1/ («) is the Whittaker 
function. Because the qualitative behavior of Wo(?) is of dumping oscillation as a function of z, the normalized 
Hawking spectrum (|47]) of present case is expected to be oscillatory with decreasing amplitude as a function of ui. To 
estimate the high frequency amplitude, note that the asymptotic form of Erfc(z) is given using that of Wp,i/(z) as 
Erfc(z) - e' z /2z, as \z\ — > 00. Consequently we find AE(0,w) - 0(e~ ul ^ cn y/cti/u) as w — > 00, then from (|53|), 
the high frequency amplitude of the Hawking spectrum can be roughly estimated to be IiiiVii Vfl w ) ~ 0(e~ 2w /w) as 
w — > 00. 



4.3.3 Case: T(r?) - R* (77) cx 77 s 

The third case is that the scale factor satisfies, 

T{rj) - R*{ij) = -ctj 3 + brj + p 



const > 



b.p = const. 



(54) 



Here note that even if there is a term proportional to rj 2 , it can be vanished by shifting the origin of 77. The differential 
of this equation gives, 



^) = ^AEM a(v) 

drj r 



-3ct7 2 



ro 



1 - 



Ra 



a(v) r o 



(55) 



From this equation we find that, if a > 0, the behavior of a(rf) is the exponential grow for 77 > yjb/ic. That is, it can 
be said that this case also corresponds to the asymptotic form of the open FRW universe at a remote future, a(t) cx t. 
As for the case of the previous subsection 4.3.2, the difference from the thermal case comes from the detailed behavior 



of a(t). For example, because c > 0, the velocity of expansion da/dr) of the present case ( |55| ) should be larger than that 
of the thermal cas e ([l9| ) for r\ > 0. 

The spectrum (|47|) of the present case is expected not to be the thermal one. To evaluate the normalized junction 
coefficient (fl6|), we should carry out the integral, J^ f drj exp[ic5l77 3 + i(u> — b£l)rj\. But it is difficult at least for the 
author to evaluate this integral. Therefore in order to guess roughly the behavior of the Hawking radiation, we try using 
the non-normalized representation (Khfl). The non- normalized junction coefficient (p0[) is exactly calculated to be, 



/((l) —iuTQ—ipQ. 



J- 1/3 



yJ\LJ -bfl\ 
[3(cn)V3]3/2 

— 6 J^i ^ 3/2 



k 3(cn)vs 

where the integral representation of the Bessel function is used 

2TT,fik 



sgn(w - bil) J 1/3 



\uj-bfl\ 



3/2' 



dxe i( - qx) +lkx 



(3g) 3 / 2 

Then the deviation factor becomes, 



J -1 



/3 



3/2' 



sgn(fc) J 



1/3 



3q 



3/2' 



for q > . 



27c n 2 



\w-bSl\ A+(n, uj) 2 + e 27ra/K (uj + b Q) A_ (fi, uj) 2 
+ 2 e nU / K ^/\u) 2 - c 2 fl 2 \ cos (2p fi) A+ (O, w) A_ (fi, w) 
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where 



J '-1/3 
J -1/3 



601 



3/2' 



3(cO) 1 /3 
cj + 60 N 3/2 * 



3(cO) 1 /3 



sgn(w - 60) J1/3 
[ / i^ + 60 

Jl/3 



601 



3(cn)V3 



3/2 



3/2' 



This result denotes that the Hawking spectrum is oscillatory as a function of u which is the frequency observed by a 
comoving observe at a remote future r\ — > +oo. The amplitude can be roughly estimated to be w (0| Nf,u |0)w ~ ujl ^ 2 
as a; — * oo, where the asymptotic form of the Bessel function is used, J v (z) ~ y/2/nz cos[z — (2i/ + l)7r/4] as |z| — + oo. 
Because this result is not normalized, the high frequency spectrum diverges. But the qualitative behavior that the 
spectrum is oscillatory may hold even for the normalized one. 



4.3.4 de-Sitter case: a(rj) = —c/rj 

Here we turn our attention from the open universe to the other universe. Consider the de-Sitter case, 

— Q 

a ( r l) = — j f° r ? / < an d c = const. > 0. 
V 

The relation, dt = a(r))drj, gives a(t) = exp(i/c), that is, this is the inflationary universe. 
From ( |lO| ) and the definition of R*(rj), we obtain the phase factor, 



I - kr 2 -R g ) In ( - 



Ra 



V cr 



P, 



where p is the integral constant due to (|10|). Hereafter in this subsection, we assume that the junction surface £ of the 
SCU is much larger than the black hole, R g <C R(r)) = —cro/rj. Then the approximation, l/rj+ R g /cr$ ~ 1/r), becomes 
valid, and the normalized junction coefficient (Eg) is expressed as 



.4 



(i) 



-iujTq —ip Q 

Vp 

-iujTq —ip O 










/ 













d?7 e 



cro O 



i(o)77-cr ii/v) rl in ( c y/ 1 - kr o 2 - R s 



due z/Hu-l/u) u -l-v _ 



where z = i 2\/cro O To, v = — 1 — iCl ( c -y/l — k Tq — R g ^j, and the variable of integration is changed as r\ — > it = 

rj yjujjcr^l. Here let us take the approximations, 77/ — ► and 77.^ — > —00, and further consider the modification, 
oj — > lim e _>o w e _le . Then the integral path can be analytically connected to the complex u-plane and deformed to C 
shown at the figure |^, which gives 



.4 



(i) 



»7p 



(56) 



where the integral representation of the Hankel function is used, iirW v (z) = — Jpdue 2 ' 2 '" -1 '"' u~ 1 ~ u for 5ft(z) > 0. 
This is valid only for uj > 0. But we can obtain the same result (56) for uj < by routine calculations with taking the 
branch cut of In it in the upper half it-plane (— 1 = e~ 17r ). Therefore the deviation factor is given as 



D H (fi,u) 



7r 2 cro 



H« (*2^ 



= 27rfi/K 



Hi 1 ) (2^ 



+2e^ K U ^ e W2-«pn 1 " H,W (t2^0^) H« (2^0^ 

This is essentially the dumping oscillation about uj. 

The amplitude of the high frequency spectrum is estimated using the asymptotic form of the Hankel function, 

Ul 1] (z) ~ v /27^Ie*( 2 -( 2l '+ 1 W 4 ), to be, H^ (i2Vcr Ow) ~ O ((Ow)- 1 / 4 x e 

Djj is the second term, and we find that (— 00, 0; uj) ~ 0(l/y/u) as cj — ► 00. 



/Qcj — o 



Then the dominant term of 
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4.3.5 Case: close to the thermal spectrum 



Finally let us examine the case which is close to the thermal case. According to the subsection 4.3.1, the Hawking 
spectrum (Ejj) should be reduced to the thermal one by a manipulation which makes the scale factor come to satisfy 
(fi"9|). For example, the slow cosmological expansion limit, a/a — > 0, would reproduce the thermal spectrum. Therefore 
we consider the case, 

T(rj)- R*{ri) = erf + bn+p as c — > 0, 

where I is an arbitrary natural number, and it is assumed that the scale factor is analytic as a function of 77 and 
consequently the phase factor T — R* is also analytic. This manipulation, c — > 0, should include the limit a/a — ► 0. 
As c — > 0, the junction coefficient is calculated as, 



A 0) 



-iujrl — ip H 

-iujry—ip n 



'If 



drj (l-tcfl^ 1 ) e ^ u - ba ^ 



UJ 




Q 





Vi 



Vp hi 



Here by treating the Kronecker du,b(i as a distribution and referring to the property of the delta function, dS(x)/dx 
—S(x)d/dx, the integral in the second term becomes, 



'if 



which gives 



.4 



(i) 



drjrj e 



I „i(u>-bn)r) 



<.bSl 



where u) — 2im/r]p, Q = 2irN/T p . This means Aq 

D H (n,Lo) = 



^ = for uj > and Q > 0, therefore we obtain 
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(i) 



x 2 



as c — > 
as c — > 0. 



In the limit c — > 0, the second and third terms of the deviation factor Dh do not contribute to the Hawking spectrum 
lH{T]i,rif\uj), and the first term of Dh is reduced to that of the thermal case. However in the non-linear case, c ^ so 
small, the second and third terms of Dh come to contribute to lH(T]i,T)f,ui). That is, we can expect naively that the 
intensity of the non-thermal Hawking radiation is stronger than that of the thermal one due to the exponential factor 
gi-ftyK m Dh, where the mass of the black hole is the same for both of the cases. 



5 Summary and discussion 
5.1 Summary 

How does a dynamical boundary condition like an expanding universe change the well-known properties of asymptotically 
flat black hole spacetimes? This paper is the first trial to treat this question, and deals with the effect of cosmological 
expansion on the Hawking radiation. The Swiss cheese universe (SCU) has been adopted as a concrete model of 
background spacetime. With using the SCU we could avoid the issue how to define the black hole in an expanding 
universe, and the mass of the black hole could be defined by the Kumar mass computed on the spatial section of the 
junction surface S of the SCU. Further for simplicity, the SCU has been assumed to be two dimensional in order to 
neglect the curvature scattering of a massless scalar field of minimal coupling. It has also been assumed that the SCU 
was eternal, that is, no big bang and big crunch singularities would appear and the scale factor a(rj) could be defined in 
the infinite interval of the conformal time of the FRW-side, —00 < 77 < 00. 

We have introduced a massless scalar field $ of minimal coupling, and ignored the back reaction to the background 
SCU. The junction condition of $ has required no potential confined on S. Therefore $ was not reflected at E. 
Consequently, since we have neglected the curvature scattering, a mode function which is a monochromatic out-going 
one in the BH-side has been connected at E to a mode function in the FRW-side which is of a superposition of the 
out-going modes without including any in-going one, and so has been the in-going modes. 
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The scalar field has been quantized with the WX-mode, the TR-mode and the Cfl-mode. Because the event horizon 
bifurcated in the eternal SCU, we could make use of the discussion given for the asymptotically flat eternal Schwarzschild 
spacetime, and obtained the Bogoljubov transformation from the WX-mode to the TR-mode Q. For finding the 
Bogoljubov transformation from the TR-mode to the Cfl-mode, it has been essential that the junction coefficients of 
the mode functions determined at £ could be re-interpreted as the Bogoljubov coefficients. Finally by the successive 
Bogoljubov transformation from the WX-mode to the Cfl-mode via the TR-mode, we have obtained the Hawking 
radiation (|dj) as the particle of the Cfl-mode created on the vacuum state of the WX-mode. This Hawking radiation 
was measured by a comoving observer at a remote future rj — > +00. 

However, because the back reaction was ignored, the Hawking spectrum ([ll]) has diverged as mentioned at the 
beginning of the section 4.2. This divergence arose from the infinite time interval of the observation. Therefore we have 
introduced the finite time interval normalization which restricted the observation time, T]i < r\ < r)f, and obtained the 



normalized Hawking spectrum (47) 



5.2 Discussions 

It is easily expected that the Bogoljubov transformation from the TR-mode to the Cfl-mode causes the difference of the 
resultant spectrum from the thermal one. Note that this Bogoljubov transformation does not include the cosmological 
particle creation, since the massless scalar field of minimal coupling in two dimensional spacetime is not scattered by 
the background curvature. Therefore the deviation factor Dh of the resultant spectrum ( f42| ) includes only the effect 
of the relative motion of the comoving observer to the TR-system. One may think that the resultant spectrum would 
be just a thermal one which would receive a red shift to a lower temperature due to the comoving motion. However 
it is obvious from (^) that the resultant Hawking spectrum is totally different from a thermal one, except the case of 
subsection 



4.3.1 



The reason why the spectrum is not generally a thermal one, is in the junction condition of at S. 
If a mode function can be a monochromatic in both of the BH-side and the FRW-side, we can understand the junction 
condition as a simple Doppler effect. However as mentioned in the previous subsection, a monochromatic out- going 
mode in the BH-side goes through S into the FRW-side, then it turns to a superposition of out-going modes in the 
FRW-side. Therefore a monochromatic out-going mode received by a comoving observer in the FRW-side, should be 
a superposition of out-going modes in the BH-side when it would be traced back into the BH-side. This denotes that 
various frequencies Q of the TR-mode are included in each monochromatic mode of frequency to in the FRW-side. That 
is, so many Doppler effects take place at the same time in one Cfl-mode, which is represented by Dh- We interpret such 
the superposition of Doppler ef fects a s the pure effect of the acceleration due to the cosmological expansion. 

According to the subsection 4.3.1, it seems that the thermal Hawking spectrum is of a special case where the second 
and third terms in the representation ( |42] ) of Dh vanish, and that the Hawking spectrum (^?j) is generally different from 
the thermal one. However as examined in the subsection 4.3.5, a non-thermal spectrum is reduced to the thermal one by 
an appropriate manipulation, for example, by the slow cosmological expansion limit, a/a — > 0. About the non-thermal 
Hawking spectrum, the subsections 4.3.2, 4.3.3 and 4.3.4 give the conclusion that the qualitative behavior of the non- 
thermal spectrum is the dumping oscillation as a function of the frequency u> measured by a comoving observer. Further 
as discussed at the end of the subsection 4.3.5, the intensity of the non-thermal Hawking radiation is stronger than that 
of the thermal one. Then we can describe a picture from the viewpoint of the black hole thermodynamics that a black 
hole with an asymptotically flat boundary condition stays in a lowest energy thermal equilibrium state. When a black 
hole is put into a dynamical boundary condition, it is excited to a non- equilibrium state, and emits its mass energy with 
stronger intensity than the thermal one. However, about our general aim given at the beginning of section [l], we could 
not search for a concrete example of the non-gravitational and non-equilibrium system which corresponds to a black 
hole in an expanding universe. 

So far we have assumed two dimensional SCU. In extending to four dimension, there arise three problems due to 
the curvature scattering: (i) the initial vacuum state, (ii) the junction of $ at £ and (iii) the cosmological particle 
creation. For the first issue, as mentioned in the subsection 2.2.3, we should modify the vacuum state on which the 



particle creation is estimated. The second issue means that, even if F = at (|69| ), a reflection of $ due to the curvature 
scattering takes place at E, that is, a monochromatic out-going mode in the BH-side should be a superposition of the 
out-going and in-going mode s in the FRW-side. Both of these issues (i) and (ii) do not cause any change for our strategy 
explained in the section However about the third issue (iii), we should add another step to our strategy. The new 
step is an extra Bogoljubov transformation from the Cfl-mode of initial time to that of the final time, which raises the 
cosmological particle creation. But such a Bogoljubov transformation will take a complicated form, and we can hardly 
estimate its form since it is difficult to solve the mode functions with considering the curvature scattering. However if 
one wants to research the effects of only the cosmological particle creation, it is the easiest model to consider the massive 
scalar field on two dimensional SCU [|10|. 



Finally let us discuss about the assumption of the "eternal" SCU. As mentioned at the end of section 2.1, though the 
evolution of the scale factor is mathematically given by the equation (|6^) , we dared to assume that the scale factor was 
a positive definite function a(r\) > for — 00 < 77 < 00. With this assumption we could obtain the form of the junction 
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coefficients A^ 1,2 ^ in the simple representation of (f23|). Then what should we modify if the equation ( p6| ) is taken into 
account? Because of the initial singularity a(t = 0) = 0, we should solve the junction conditions of the scalar field <f> for 
the interval i]q < r/ < oo, where r/o is given by a(r]o) tq = R g . This modifies the completeness (|lj) to 



dv r F , u (V) If,w = tt * (a; - w') + ie« u -">*> pv 



'/(> 



where pv is the Cauchy's principle value. The second term may raise some complicated term in the equation (p3|), which 
would make the following computations intricate. However if we can assume that 



dw ' pv i^L^,2) =0j (57) 



U) — W 



then we can obtain the junction coefficients of the form (|23|) with modifying the interval of integration of I a and J a 
from — oo <?7<ooto?7o<?7< oo. Further even if this modified junction coefficient makes the Hawking spectrum 



diverge, we can expect well that the same normalization method discussed in section 4.2 holds. Here, by considering the 
four dimensional case, we can find the validity of the assumption (p7[). In four dimensional case, we can also separate 
the variables on the FRW-side, <& = /f(t?) hp(r*)Yi m (9,p), where Yi m is the spherical harmonics. It can be easily 
found that the portion /f(?7) satisfies a Sturm-Liouville type ordinary differential equation. This implies that there is a 
complete orthogonal set {/f,w} satisfying ( jilt ) with appropriate range and measure about r\ and uj. Therefore it is not 
so bad to assume the background spacetime is the eternal SCU in order to extract the essence of the Hawking radiation 
in an expanding universe. 
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A Four dimensional swiss cheese universe 



The swiss cheese universe (SCU) || is the spacetime including a spherically symmetric black hole in an expanding 
universe, and obtained by connecting a Schwarzschild spacetime with a dust-dominated Friedmann-Robertson- Walker 
(FRW) spacetime at a given spherically symmetric timelike hypersurface, E, by the Israel junction condition with no 
energy density confined on the junction surface ||. In this appendix, we show a sketch of constructing the SCU in four 
dimension. Hereafter as the terminology, let the word "BH-side" denote the spacetime region inside E where the metric 
is given by ([58]), and "FRW-side" for the region outside E where the metric is (|59|). 
The metric of a Schwarzschild black hole is 



as BH 



Hi 
R 



dT 2 



R 



dR 2 + R 2 (d0 2 + sin 2 6dip 2 ) , 



(58) 



(T, R, 9, ip) is the Schwarzschild coordinate and R g = 2GM, M is the mass of the black hole and G is the 



where 

gravitational constant. The metric of a FRW universe is 



ds% 



-dt 2 + a(t) 2 



dr 2 



1 — kr 2 



de 2 



sin 9dip 



(59) 



where x F = (t,r,9,ip) is the comoving coordinate, a(t) is the scale factor and k = ±1, is the spatial curvature. We 
can make the angular coordinates 6 and tp to be common to both of these coordinates due to the spherical symmetry. 
Further we assume that the Friedmann equation, that is the evolution of a(t), is of dust-dominated one, 



k 8irG po 
_ _ 



a" 



(60) 



where po (= const.) is of free parameter at this moment, but will be fixed consistently in considering the junction 
condition later. 

In constructing the SCU we choose E to be spherically symmetric and timelike. Therefore the angular coordinate 
9 and p can be chosen as the spatial coordinates on E, and it is possible to use t as the temporal coordinate on E. 
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That is, the coordinate we use in E is expressed as u % — (i, 6, <p). [] The radial coordinates of E measured from BH-side 
and FRW-side are generally given as a function of t, like R(t) and r{t) respectively. However especially about r(t), we 
can specify it to be constant, r = tq = constant. To understand this feature of radial coordinate of S, it is important 
to notice that the mass of the black hole M is assumed to be constant. If r(t) 7^ constant, there should be an energy 
flow of dust-matter through E, since the dust-matter rests on the comoving coordinate. Consequently, because it is also 
assumed that no energy density is confined on E, the source of the energy flow should be the mass of black hole, that is, 
M should have time dependence. [] This contradicts M — constant. Hence r(t) = r = constant. On the other hand, 
another function R(t) cannot be specified at this moment, but will be determined by the junction condition. 
The Israel junction condition at the junction surface E with no energy density confined on E, is given by 



(61) 
(62) 



where hs_f/|s and ^bhIh are respectively the induced metric and the extrinsic curvature of E measured from the BH- 
side, and similarly to hp\s and K F |s with respect to the FRW-side. The components of h and K on E are respectively, 



Hj — miv and Kij = n^ u x r ' { Xj, wnere x r ^ = ox^/Ou- , g^ u 

normal to E. In the FRW-side, the unit normal to E is obtained from the expression of E, Sf = f 



where x 



dx^ jdu % 



is the metric of four dimension, 



n M is the unit vector 
r = 0, 



,A» 
F 



1 



v\Sf,oiSf I 

In the BH-side, Sbh = R — R(t) = gives E. Using the the junction condition (|6l|), we obtain the unit normal to E 



(63) 
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f(t) ,0,0 



(64) 



The results of the junction conditions (|61j) and (|62[), are summarized into three independent equations. One of them 
specifies the location of E in the BH-side, R(t) — a(t) ro in BH-side. The other two equations give the form of the 
function T(t) and the evolution of the scale factor on E, 



dT(t) _ s]\ - kr 2 



dt 



1 - Rg/R(t) 

k Rg 1 



Comparing the result (J66|) with the assumption (pOJ) , the consistent choice of the constant pg is 



M 



Po 



(4/3) 7 



(65) 
(66) 

(67) 



To draw a picture of the SCU, assume R(t c ) > R g holds at an arbitrary time t c for the case that the FRW-side is 
open k — — 1 or flat k = 0. Then, because a(t) increases monotonically in these cases, the relation ( |67[ ) provides us with 
the following picture of the SCU: in a dust-dominated FRW universe including no black hole, the dust-matter contained 
in a spherically symmetric region of comoving radius ro starts to collapse, and finally a spherically symmetric black hole 
of mass M = (4/3)7rrQ po is produced in an expanding universe. The conformal diagram of this interpretation is shown 
at the figure [|, which makes it clear that a black hole can be defined in the SCU with open and flat FRW-side. On 
the other hand, it is impossible to define a black hole for the case of closed FRW-side. Because a(t) decreases to zero 
in a finite future, no causal curve cannot avoid encountering the big crunch singularity, hence every spacetime point is 
included in the causal past of the singularity and we cannot define a black hole. 

Finally as the definition of the mass of the black hole in the SCU, we adopt the Komar mass calculated on the 
junction surface E. The general form of the Komar mass Qk @ is 



Qk = 



8nG 



t)A 



where A is a spacelike region, dA is the boundary of A, £ is the timelike Killing vector in A and dS^ is the surface 
area form on dA. Using the Einstein equation and the Killing equation, it can be proved that Qk is invariant under 

5 It is also possible to use T as the temporal coordinate on S. The reason we chose t is just the later convenience. 

6 For more correct discussion, consider the Komar integral evaluated on an arbitrary closed twa dimensional spatial surface including the 

Satial section of the event horizon. Then the resultant Komar mass is exactly M for the metric 1581) as calculated at the end of this appendix 
This means that the energy flow through S affects directly the mass of the black hole, if S is located outside the event horizon. 
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Figure 4: Conformal diagram of the SCU interpreted as the collapse of the dust-matter in flat or open FRW universe. 

any deformation of A if £ lasts to exist on the deformed region. For the SCU, dA is the spatial section of E at T = 
constant, and we can let £ = <9t, then the Komar mass is obtained to be Qk = M. Here note that an important issue 
is left: the normalization of our Killing vector £ = dx is not given at the junction surface nor at the spatial infinity of 
the FRW-side, but given at the asymptotically flat region of the full Schwarzschild spacetime which is not present in 
the SCU. Therefore we should say that, with accepting the vector dx as a timelike Killing vector in the BH-side, the M 
appeared in ( |58| ) can be understood as a mass of the black hole in the SCU which is invariant under any deformation of 
the spatial hypersurface in the SCU. Further as a by-product of this normalization, the surface gravity k of the event 
horizon is evaluated to be k = l/2R g . 



B Junction condition for scalar field 

The junction condition for connecting a scalar field <p at a given hypersurface £ is discussed in this appendix. As a 
terminology, we call the spacetime regions separated by £ A-region and B-region, and denote a quantity Q evaluated 
in A-region and B-region as Qa and Qb respectively. We consider the four dimensional background spacetime in the 
Gaussian normal coordinate with respect to a hypersurface E, and the scalar field 4> with an arbitrary potential V{4>), 
The metric is expressed as 

ds 2 = dx n 2 + Qij dx l dx^ , 

where x % is the coordinate intrinsic to the hypersurface E, x n is the coordinate vertical to E, and — g nl1 is the unit 
normal vector to E, that is, n 2 = 1 if E is timelike or n 2 = — 1 if E is spacelike. Hereafter we set the surface E is placed 
at x n = 0. The Klein- Gordon equation is 

By integrating this equation along x n direction through E, we find 

dx n \3<p= [ dx n — . 

, J-e d<t> 

On the other hand by definition of the d'Alembertian, we obtain 

J dx n D(t>= J dx n d n d n cl>+ J dx n (r^d^ + V.V^) — ► d n $ A -d n (\> B ase^O, 

where d n = n^d^cf), and it is assumed that the A-region is the region of x n > and B-region is of x n < 0. Here assume 
that there is a potential confined on the surface E such as 

^ -"<*•>■ 

where F denotes a "surface potential force" derived from the potential confined on E. Further it is natural to require 
the continuity of <f) on E. Hence finally we obtain the junction condition of a scalar field on a given surface E to be 

0a|b = ^b|e (68) 
n'Xd^Als - n^d^Bls = F, (69) 

If no potential is confined on the surface E, the surface force vanishes F = 0. This junction condition of a scalar field 
is very similar to the Israel junction condition for connecting geometries of spacetimes in the sense that a kind of 
source of the field under consideration is required in order to satisfy the field equation even on the junction surface. 
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C Second term of equation (|27|) 

In this appendix we discuss how the second term of ( p7j ) vanishes and what it means. To begin with, we need to know 
the behavior of the phase —itl[T(r)) — R*(rj)]. Because the background is the eternal SCU, the junction surface £ 
reaches the past temporal infinity as rj — » — oo and the future temporal infinity as r\ — > +00. Then the equations (|3|) 
give T(rj) — R*(rj) — » ±00 as r\ — > ±00. 

Note that the mode functions V'bh si should be the basis of the general solution of the Klein-Gordon equation, 

/oo 
-00 

where the coefficients depend only on f2. This means that we should consider the property of the modes V'bjj n 

under the integration with a function of fl. Let iJ(fi) be a test function, we can calculate as f\ 

dfl c -m\T(n)-R*M] 



H(Q) 



-in[T( v )-n*( v )] 



T( V )-R*(r,) 
then because T(rj) — R*(rf) — > ±00 as 77 — > ±00, 
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i I d H(Q) 

t(v) - R*(v) [dn 



-in[T( v )-R*( v )] 



an- 



-in[T {v )-R*{ v )] = Q asr] ^ ±OQ _ 



Hence the second term of (^) vanishes ^ whenever the mode functions V'bjj n are considered under the integration with 
respect to f2, 



sgn(?7) 7T 
(4 7 r) 3 / 2 v /|fi| 



r)—-\-oo 



= 0. 



(70) 



Further from this convergence, the equations ( p6| ) denote Ia(£1,u) — —Ja(Q,(j) whenever the mode functions V-'bh si 
are considered under the integration with respect to Q. 



[3 

[5 
[G 

ir 

[8 
[9 
[10 



References 

[1] Barrow J.D., 1992 Phys.Rev.T)46 R3227; Barrow J.D. and Carr B.J., 1996 Phys.Rev.D54 3920; Shibata M. and 
Sasaki M., 1999 Phys.Rev.B60 084002; Jacobson T., 1999 Phys.Rev.Lett.83 26 99; Saida H. and Soda J., 2000 
Class. Quant. Grav.17 4967; Harada T., Carr B.J. and Goymer C.A., 2001 Preprint |astro-ph/0102058| (submitted to 
Phys.Rev.B); Sakai N. and Barrow J.D., 2001 Class. Quant. Grav. 18 4717 

[2] See for example and the references therein, Wald R.M., 1994 "Quantum Field Theory in Curved Spacetime and 
Black Hole Thermodynamics", the University of Chicago Press 



Hawking S.W., 1975 Commun.Math.PhysAZ 199 

Hawking S.W. and Ellis G.F.R., 1973 "The large scale structure of space-time", Cambridge University Press 
Dyer C.C. and Oliwa C, 2000 Preprint ^stro-ph/0004090| 
Israel W., 1966 Nuovo Cim.44B 1, 1967 Nuovo Cim.48B 463 



See for example the chapter 5 in the good review article: Townsend P.K. 1997 Preprint gr-qc/9707012 
Hartle J.B. and Hawking S.W., 1976 Phys.Rev.D13 2188 
Unruh W.G., 1976 Phys.Rev.D14 870 

See for example the section 3.4: Birrell N.D. and Davies P.C.W., 1982 "Quantum fields in curved space", Cambridge 
University Press 



7 Mathematically this is the Riemann-Lebesgue Integral Theorem. 
8 In an exact word, this is the weak convergence. 



21 



